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ABSTRACT 

In a line caustic crossing microlensing event, the caustic line moving across the surface 
of the source star provides a direct method to measure the integrated luminosity profile of 
the star. Combined with the enormous brightening at the caustic crossings, microlensing 
offers a promising tool for studying stellar luminosity profiles. We derive the amplification 
behavior of the two extra images that become partial images conjoined across the critical curve 
at a line caustic crossing. We identify the multiplicative factors that depend on the caustic 
crossing point and the relative size of the star, and the shape function that depends on the 
stellar luminosity profile. We examine the analytic limb-darkening models - linear, square 
root, and square - using the analytic form of the shape function. We find that the microlensing 
lightcurves must be determined to an accuracy of better than 0.3-0.8% in order to be able to 
determine the linear limb-darkening parameter c\ with precision of 5c\ = 0.1. This is similar 
to the accuracy level required in eclipsing binaries as reported by Popper (1984). 



Subject headings: gravitational lensing - binary stars 
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The caustic curve of a binary lens is where the number of images changes by two, and the two images 
that appear or disappear are extremely bright. The caustic curve is a geometric curve with zero width, and 
the differential limb darkening profile of a lensed star (of a typical size of ~ 1 — 10/uas) can be readily 
scanned during a caustic crossing. The achromaticity of gravitational lensing and the large magnification of 
the two extra images help microlensing to be a rather clean tool for studying the luminosity profiles of the 
lensed stars. 

The microlensing method has several advantages: 1) The universality of gravitation guarantees 
that any star can be lensed, and its limb darkening can be studied through microlensing . That is not 
always the case in other methods. For example, the lensed star of the recent line-caustic crossing binary 



microlensing event MACHO-98-SMC-1 (|Rhie et al. (1999); Afonso et al. (2000) and references therein) 



was a (spectroscopically) normal low metalicity A star ( Albrowetal. 1999| ). Normal A stars have not been 



found in short-period binary systems ( |Abt 1993[ ), and their limb darkening is unlikely to be studied through 
eclipsing. 2) The auxiliary object that accommodates high resolution scanning of the lensed star is a 
foreground (lensing) object that has no physical interaction with the lensed star under study. Combined with 
enormous brightening, the physical independence simplifies the interpretation of the event. In an eclipsing 
binary, the very close interaction of the eclipsed star with the companion (the auxiliary object) requires an 



elaborate modeling of the binary system with many parameters ( popper 1984| ) and makes it more difficult 
for the method to be robust. 3) Microlensed stars are in the Galactic bulge (~ 8 kpc) or in the Magellanic 
Clouds (~ 50, 60 kpc), and microlensing allows studies of the luminosity profiles of faint distant stars. 

Resolving stars without auxiliary means has not been feasible technologically except for the biggest 
stars in the s ky: the Sun and more recently Betelgeuse (M2 supergiant) studied with HST (G illiland & 
Dupree 1996). Most of the studies of stellar luminosity profiles have been done with a handful of detached 
eclipsing binaries. In eclipsing binaries, the companion traverses in front of the star and effectively 
scans the luminosity profile of the eclipsed star by absorbing the light. How well one can determine the 
profiles depends on the well depth of the eclipse light curve, and it has been estimated through numerical 
simulations that 100 or more observations within minima with a relatively small photometric error of 
~ 0.5% are necessary to be able to determine the linear limb darkening parameter of the eclipsed star with 



the uncertainty of 5c\ = 0.1 or less ( [Popper 1984J ). These simulations were done with simplified models, 



and obtaining the precision of 5c\ =0.1 may turn out to be more demanding in practice. 

In a caustic crossing gravitational microlensing, the caustic curve traverses across the surface of the 
star and scans the integrated luminosity profiles of the star by partially magnifying the star. This partial 
magnification is due to the two extra images, and the magnification behavior of them can be expressed 
as a simple analytic function integrated over the luminosity profile. We examine analytic models of the 
luminosity profiles and find that relatively small photometric errors are required for the determination of 
the limb-darkening profiles: In order to be able to determine the linear limb darkening parameter c\ with 
an uncertainty of 8c\ = 0.1 or less, photometric accuracy of 0.3 — 0.8 % is required depending on the value 
of ci (see figure [7] and |8]). 

It has been a common practice to incorporate the effect of the luminosity profile of the lensed star for 
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the high resolution behavior of the light curve (for example, supernovae lensing by Schneider and Wagoner 



(1987)). For stable stars, the linear limb-darkening law (Milne 1921; Gray 1995) has been the default model 
for stellar luminosity profiles. 

1(h) 

——=c + cih; c +ci = l, (1) 
where h is the cosine of the angle between the line of sight and the normal vector of the surface element 



of the photosphere. The Sun is known to follow the linear law in the optical bands (Gray 1995). We had 



adopted c\ = 0.6 as a generic value for main sequence stars in the calculation of the detection probability 



of microlensing earth mass planets (Bennett & Rhie 1996). However, the linear models are suspected to 



be inadequate especially for hot stars, and non-linear terms have been added in the functional form of 
the analytic models. Recently, Diaz-Cordoves, Claret and Gimenez studied LTE (local thermodynamic 
equilibrium) model atmospheres (with standard solar compositions) calculated by Kurucz (1991) in 



search for suitable analy tic limb-darkening models (piaz-Cordoves & Gimenez 1992& D iaz-Cordoves, 



Claret, & Gemenez 1995; |Claret, Diaz-Cordoves, & Gemenez 1993 ) that fit the numerical results. The 



limb-darkening models in their work can be written in the following form. 

1(h) 



lo 



Y,cnh n - ^c n = l: n = 0,0.5,1,2 (2) 



A homogeneous model is given by c n = 0; ti/0,a linear model is given by c n = 0; n / 0, 1, and a 
non-linear model refers to a combination of the linear terms and one of the non-linear terms (c2 / 0: square 
model; C0.5 7^ 0: square-root model). The larger the index n, the more darkened in the limb is the radial 
luminosity profile given by h n . I is the center luminosity (at h = 1), and the condition J2 c n = 1 ensures 
that the linear combination on the RHS of the equation (|2|) has value 1 at h = 1. 

In a binary lens, a source star outside the caustic curve produces three normal images and a source 
star inside the caustic curve produces extra two bright images in addition to the three normal images. At a 
caustic crossing, the surface of the star is divided into two regions by the caustic line moving at a constant 
velocity, and only the part of the star that falls inside the caustic produces two extra images. These partial 
images are connected at the critical curve (of the lens equation), and this criticality is what underwrites 
microlensing as a potential tool for measuring stellar luminosity profiles. 1) The magnification of these 
two images is governed by the generic critical behavior of lensing, and this generic behavior simplifies the 
interpretation of the lensing. Thus, the light curve of the two partial images bears direct measurements of 
the integrated luminosity profile of the source star (measured in time through the moving caustic line). 
2) The critical curve is the loci of the images of divergent magnification of point sources. Thus, the two 
partial images joined across the critical curve are highly magnified and dominate the behavior of the total 
light curve. The three normal images are full images, insensitive to the luminosity profile of the source star, 
and their brightness varies relatively slowly in time (unless the crossing is very near a cusp). 

The nomenclature critical curve derives from the fact that the lens equation is critical, or, stationary on 
the curve. One of the eigenvalues of the Jacobian matrix of the lens equation vanishes on the critical curve, 
and the lens equation is quadratic in the eigendirection of the vanishing eigenvalue. As a consequence, 
the amplification of the images in the neighborhood of the critical curve is proportional to the inverse 
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square root of the distance of the source to the caustic. The proportionality constant that determines the 
overall peak amplification of the light curve of the two extra images is given by the inverse square root 
of the derivative of the Jacobian determinant in the critical direction (the eigendirection of the vanishing 
eigenvalue). General discussions of this critical behavior of lensing can be found in the literature (e.g., 



Schneider, Elhers, & Falco (1992)). However, we include a brief derivation of the critical behavior of the 
class of binary lenses not only because the class is directly relevant to experiments, but also because its 
simplicity allows easy access to concrete understanding of the critical behavior. 

If the complex variables z and w denote an image position and its source position respectively, the 



binary lens equation is given as follows dBourassa, Kanotowski, & Norton 1973| ; pennett & Rhie 1996 
Rhie et al. 200C| ). 



oj = z-R 2 E (^— + ^-) ■ ei + e 2 = l (3) 

\Z — X\ z — x 2 J 

where e\ and e 2 are the fractional masses of the lenses located at x\ and x 2 respectively. The lens equation 
is intrinsically a relation between angular position variables, and the lens plane (the complex plane 
parameterized by z or u) can be put anywhere along the line of sight. The lens plane is a linear plane normal 
to the line of sight and is equivalent to the angular space because scattering (deflection) angles involved in 
microlensing are small (<< 1). If the lens plane is considered to pass through the center of the mass of 
the lenses, the Einstein ring radius of the total mass M of the binary lenses is given by R E = y/AGM D 
where D is the reduced distance: D^ 1 = D^ 1 + D^ 2 , and D% and D 2 are the distances along the line of 
sight (or in the radial direction) of the observer and the source star from the lensing system respectively. 
(It is conventional that "Re" is reserved for the Einstein ring radius given by y4GMD.) The position 
variables can be rescaled by Re so that the unit distance scale of the lens plane is set by the Einstein ring 
radius: Re = 1. (In the following, all the position variables are measured in units of Re- However, 
we may replace 1 by Re when we find it useful for clarity.) The so-called image plane is the lens plane 
parameterized by the image position variable (here z) and the source plane is the lens plane parameterized 
by the source position variable (here u>). 

The linear (differential) behavior of the lens equation is given by 

(4) 






where Zi = z — X{ : i = 1,2, and 

h 

The eigenvalues of the Jacobian matrix J are 



= Ule**. ( 5 ) 



A± = 1 ± \k\ , (6) 

and the corresponding eigenvectors are 
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Where \k\ = 1, the eigenvalue A_ vanishes, and the Jacobian matrix J behaves as a projection operator. 
Thus, the curve defined by \k\ = 1 is called the critical curve. Obviously, the Jacobian determinant 
J = 1 — \ k\ 2 vanishes on the critical curve. If dz is a displacement from a point z on the critical curve, duo is 
a non- vanishing displacement vector parallel to the upper component of the eigenvector e + with eigenvalue 
A + = 2. Let E+ and E- denote the upper components of the eigenvectors (the lower components are their 
complex conjugates), i.e., E + = e~ tLp and E- = ie~ lip . Then, dz = dz + E + + dz-E- where the linear 
coefficients dz± are real, and, 

duo = \ + dz + E + = 2 dz + E + = 5uo\ + E + . (8) 

Because of this projection effect, the corresponding area element (in the w-plane) vanishes, i/2 duoAduo = 0, 
and the magnification of its image diverges: | J| _1 = dz A dz/duo A duo = oo. Now, the caustic curve is, by 
definition, a curve on the w-plane which the critical curve is mapped onto under the lens equation. Thus, 
if dz is an arbitrary displacement tangent to the critical curve, the induced displacement duo in the o;-plane 
is tangent to the caustic curve. Since the induced displacement duo is always in the direction of E + as we 
have seen in equation ([|), E+ is tangent to the caustic curve everywhere; or, the critical direction E- is 
normal to the caustic curve. Therefore, when the source star approaches the caustic curve, the distance to 
the caustic S can not be expressed as a linear function of the image displacement from the critical curve dz. 
If 5uo2 is the second order contribution, 

8uo 2 = \{~dJ) {dz+ - i dzJf . (9) 

If 8uo2 = 5uo2+E + + 5uo2-E-, 

Su 2 - = - (d+J (2dz+dz-) - d-J(dz\ - dz 2 _)} , (10) 
where d + J and <9_ J are the derivatives of the Jacobian determinant J in the eigendirections. 

dJ = dz + d + J + dz^d-J . (11) 

Since J = on the critical curve, the Jacobian determinant near the critical curve is simply J = d J in the 
linear approximation. Now, we need to solve the lens equation to find dz + and dz- for a given displacement 
5uo. If we consider only the leading order terms, 

5uo + = 5uoi + , 5uo- = 5uo 2 - , (12) 

then, dz + = 5uo + /2, and dz_ is the solution of the following quadratic equation. 

duo- = ^ (d-J dz 2 _ + d+J 5uo+dz_ - d- J^j-^j (13) 

Solutions exist where the discriminant is non-negative. 

A5uo_ d-J+ \8J\ 2 5uo\ > . (14) 
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The equality holds on the critical curve (dJ = 0) as one can verify from equations (11) and (13). The 



Jacobian determinant can be calculated for the two solutions using (11). 



J = dJ = ± v /4d_J 5u- + \dJ\ 2 Su% (15) 

Now, let's note from the equality in ( JT4| ) that the caustic curve is quadratic: Su- oc Sufi (if 9_J / 0). 
Therefore, the eigenvectors, E + and EL, rotate ("frame rotation") as they move along the caustic curve. 
In this lowest approximation, however, we can ignore the rotation, and the displacement of the point 
{5ui+ , Su>- ) from the caustic curve can be approximated by 

= 5u- - 6u-(J = 0) = 5u- + j^j—r • (16) 
The distance is 5 = | Aw_ | , and the Jacobian determinant takes a simple form. 

J{5u + ,5u-) = ±J4dU . (17) 

The two images have opposite parities, and their magnifications A± have a simple relation to the distance 
of the source to the caustic 5. 

A± = 7WT\7I ; {s|A "-' <18) 



The condition for the existence of the two extra images ( |14[ ) can be rewritten elegantly as follows. 

<9_JAu;_>0 (19) 

Since the eigenvector £L is normal to the caustic curve, if a source position with Au;_ > is in one side of 
the caustic, then a source position with Aw_ < is in the other side of the caustic. So, the inequality ( |l9| ) 
tells us that a source in one side of the caustic curve (determined by the direction vector d-J E_) generates 
two extra images, but a source in the other side of the caustic does not. Thus, a source star produces two 
partial images at a caustic crossing. In lensing where the wave property of the photon beam can be ignored, 
the caustic curve is a geometric curve with zero width, and the source star is resolved by the caustic as a 
form of two partial images that change in time. (The side of the caustic where the two extra images are 
produced is the inside of the caustic.) 

The inverse square-root dependence of the amplification A± on the distance 5 of a source position 



from the caustic curve in equation ( |18| ) doesn't hold if the second order contribution vanishes, 5u>2 = 0. 
The source displacement from the caustic 5u>2 induced from the image displacement dz (from the critical 
curve) in equation (|l4|) vanishes, when 8 J = 0. 



dJ = d(l - \k\ 2 ) = -KdR (20) 

So, there are two cases where dJ vanishes: k = or 8k = 0. The condition k = defines limit points 
where J = 1 not the critical curve where J = 0. Therefore, 5lo2 vanishes only if 9k = which is the 



condition for bifurcation of J = constant curves (Rhie & Bennett 1999). There are two cases where the 



critical curve has bifurcation points. If we let £ = \x\ — X2I t> e tne separation between the lens elements, the 
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critical curves with bifurcation points satisfy special relations between the lens parameters I and e<i (R hie 
& Bennett 1999). The critical curve has one bifurcation point on the lens axis if £ = £ + . 

£+ = (^/il+^P ; 1 < £+ < 2 (21) 

The critical curve has two bifurcation points off the lens axis if £ = £- . 

£- = (^I+^r^ 5 ~r < I- < 1 ( 22 ) 

v2 

Since the lens equation is continuous (actually differentiable) in the neighborhood of the critical curve (the 
lens equation is singular only at the lens positions where J = — oo), the caustic curve bifurcates where the 
critical curve bifurcates. At the bifurcation points (four-prong vertex) of the critical curve, the lens equation 
is stationary along the critical curve. Therefore, the caustic curve forms a (doubly degenerate) cusp where 
it bifurcates. If <9_ J = 0, the critical curve is tangent to EL , and the lens equation is stationary with respect 
to a displacement along the critical curve. Therefore, the caustic curve forms a cusp where d- J = 0. 
On the bifurcation points, both components d- J and d + J vanish because the Jacobian determinant J is 
stationary in every direction: d J = 0. Thus, if we consider only the line caustic crossings away from cusps, 
the quadratic behavior of the lens equation in the direction of E- remains valid, and the amplification of 
the two extra images across the critical curve can be discussed using the expression in equation (filf). The 
caustic curve has 8 cusps if £ > £ + , 6 cusps if < £ < £ + , and 10 cusps if £ < £-. 

The amplification coefficient |<9_ J\~ l l 2 increases toward a cusp (cL J = 0) where it diverges. Figure 
[l] shows the "central caustic" of a binary lens with the fractional minor mass e = €2 = 0.34 and the 
separation £ = \x\ — X2I = 0.65. The "central caustic" (always across the lens axis, and here contains the 
center of mass inside) has four cusps, and \/|cL J\ has four zeros at 2ip = 0, 1.16757T, 2-k, and 2.83257T 
where 2(p is the phase angle of k = e l2if (in the Riemann sheet). The four blue crosses on the caustic 
curve mark the loci of the local maxima of \/\d- J\. Since the amplification of the "two extra images" is 
proportional to it, |cL J| -1 / 2 can be considered as the "strength" of the caustic curve. (As we will discuss 
later, the actual amplification of the two images also depends on the relative size of the source star with 
respect to the Einstein ring radius, which may be considered as a scale factor.) Note that the local minima 
of \d- J| -1 / 2 are very close to the off-axis cusps. We may loosely refer to it as that the cusps on the lens 
axis are stronger than the ones off the lens axis. In the same token, we can say that the cusp at ip = (near 
the dominant mass t\ = 0.66) is stronger than the cusp at ip = it. In the case of MACHO-98-SMC-1, the 
second caustic crossing was closer to a cusp on the lens axis (2ip = 0) than the first caustic to its nearest 
cusp (2(p = 2tt), and the caustic crossing amplification peak for the second caustic crossing is bigger than 
the first (see figure 1 and 2 in Rhie et al 1999). This lens (e = 0.34 and £ = 0.65) also has two triangular 
caustic curves off the lens axis (and reflection symmetric with respect to the lens axis - real line here), 
and each has "topological charge" 1/2. The eigenvectors E± rotate by 7r around each triangular caustic 
curve (or "trioid"). The eigenvectors rotate by 2ir around the central caustic "quadroid" in figure [l] which 
has "topological charge" 1. These off-axis caustics are "weaker" than the caustics which cross the lens 
axis. For example, the maximum of y/\d-J\ on these "trioids" (closed curve that is smooth except at three 
cusps) 6.93 which is about 7 times larger than the maximum on the "central caustic" ("quadroid": closed 
curve that is smooth except at four cusps). 
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If <9_ J = 0, equations ( JTo| ) or (13) becomes linear in dz^, and one obtains a linear solution whose 
amplification depends linearly on the inverse of the distance of the source to the caustic: ^4 C usp oc 1/5. 
When the source exits the caustic through a cusp, this image crosses the critical curve. This is the third 
image that diverges at a cusp crossing. The other two image solutions can be obtained by including the third 



order terms in the expansion of 5u)2- in (10). If the source is sufficiently small compared to the caustic 



curve, the probability of crossing the caustic through a cusp is small. Out of about ten caustic crossing 



binary lensing events discovered so far ( jAlcock et al 1999[ ), MACHO-97-BLG-28 was the only event where 



the crossing was through a cusp (Albrowetal. 199S). The source star of the event MACHO-97-BLG-28 



was a giant star. Most of the lensed stars in microlensing experiments are main sequence stars. Here, we 
consider caustic crossings away from cusps only. 

If the source is sufficiently small, the curvature of the caustic curve can also be ignored, and we can 
assume that the luminosity variation of the light curve over the caustic crossing is largely determined by 
the critical behavior of the lensing amplification and the luminosity profile of the source star. If we let the 
line caustic be along the y-direction, the lensing amplification is independent of y, and we only need to 
know the "1-d luminosity profile" of the source in the normal direction, say x-direction, which is obtained 
by integrating the 2-d profile along the y-direction. It is worth emphasizing that the luminosity profile we 
probe directly in a line-caustic crossing microlensing event is the 1-d luminosity profile of the source star. 
(The "irregularities" of the stars such as star spots are not addressed here. The 1-d luminosity profile will 
depend on the position of the spots with respect to the direction of the caustic line.) 

If the size of the star is r* (in units of Einstein ring radius), the luminosity profile in equation (||) 
can be rewritten in terms of the radial coordinate r of the stellar disk using h 2 = 1 — r 2 /r 2 r We let 
r 2 = r 2 (x 2 + y 2 ) such that x and y are the Cartesian coordinates of the unit stellar disk, and define 
a dimensionless 1-d luminosity function f(x) = J2 n °nfn{x) based on the radial luminosity profile in 
equation (^). 

r f\/l-x 2 . . n n+l 

f n (x) = J h n dy = J (l - r 2 /rl) 2 dy = 2b n (l - x 2 ) 2 , (23) 

where b = 1, 6 . 5 = 0.87362649, h = vr/4, and b 2 = 2/3. 

b n = [ (l-t^dt 



o 



The integrated flux is given by a byproduct of b n 's. 



F n 

We find 6 L5 = 0.71849148 and 6 3 = 3vr/16, and 

F = vr ; F .5 = 2.51077224 



J h n dxdy = J f n dx = Ab n b n+1 (24) 



4tt 2ir tt 

— ; Fi = — : Fi = — 
5 3 2 



Figure H| shows normalized radial luminosity profile h n n/ F n and 1-d luminosity profile f n ir/F n . We have 
chosen the total flux to be tt for easy recognition of the peak values of the radial luminosity profile in the 
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plot: the peak values are 1, 5/4, 3/2, and 2 for n = 0, 0.5, 1, and 2 respectively. For the 1-d luminosity 
profiles, the peak values are given by 7r/26 n+ i: 2.00, 2.19, 2.36, and 2.67. Because of the integration 
in y, the limb-darkening effect is more accentuated in the 1-d luminosity profile compared to the radial 
luminosity profile of the stellar disk. 

The amplification profile of the two extra images at a caustic crossing can be calculated using the 
luminosity profile f(x) and the analytic form of the point source amplification function A± in equation 

©• 

UA±. = G ; A = A + + A_ (25) 
J fdx 

Since the point source amplification function A vanishes in one side of the caustic line, we may specify the 
segment of the one dimensional star (as defined by the 1-d luminosity function) that is inside the caustic, 
say, as given by an interval [r*a, r*] where (a < 1). Then, A may be written as a function of x accordingly. 

A= J— — * — = A C A*A X (26) 
V |cL J\ V x - a 

The first factor A c is determined by where on the caustic curve the crossing occurs. Figure [T] shows A~ x 
for the central caustic curve of a lens with (e, £) = (0.34, 0.65). In this case, in a line caustic crossing, 
A c ~ 1-10. The other triangular caustics we mentioned above but not shown here are "weaker" and 
A c ~ 0.15-1 for a line caustic crossing. The second factor A* depends on the size of the star. 

r* = Dia* = - = — , (27) 

D\ + D 2 ctE 

where a* and i?* are angular (apparent) and the physical sizes of the star respectively, and ue = Re/Di is 
the angular size of the Einstein ring radius. Here, oe = 1/D\ because we set Re = 1. If ag ~ 1 mas, and 
a* ~ 1 - 10/xas, then A* ~ 30 - 10. 

If Cn = 1 for a given n, the amplification profile function G is determined by one function f n , and we 
may define G n = G| Cn=1 . Then, for an arbitrary luminosity profile, G is written as a linear combination of 

G n 's. 



En.C n G n (a) 



G(a) = :" w : c n = c n F n ; ^ c n = 1 (28) 

The appeai - ance of the "new" coefficients c n is due to that the directly measurable quantity in lensing is the 
amplification of the images where the normalization is the total flux while the coefficients c n 's have been 
defined based on the luminosity profile where the more relevant unit quantity is the peak luminosity. Since 
we need to be able to compare the microlensing measurements with theoretical calculations as well as 
measurements from other methods, it should be useful to consider the linear coefficients c n of the analytic 
luminosity profile functions as the basic parameters and feel free to redefine linear coefficients based on 
them as the needs arise. 

Since the amplification factors A c and A* contribute to the amplification of the two extra images as 
overall scale factors, we define g n {a) that depends only on A x (a) such that G n (a) = A c A*g n (a). Then, 
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G(a) = A c A*g(a) where the shape function g(a) is given as follows. 

9(a) = E ^" gn(a) (29) 

In an observed lightcurve, the "time" a in units of the (relative) stellar radius r # should be replaced by the 
actual time variable t. If the caustic crossing occurs at t = t cc (when the center of the star crosses the 
caustic line: a = 0), i* is the stellar radius crossing time, and i is the incident angle of the source trajectory 
(the small angle between ±E- and the source trajectory), 

a= (t - tcc) C ° SL . (30) 

So, the time axis of an observed lightcurve depends on the stretch factor 

U «* r*R E R* B\ 



cos i fi cos l v± cos i v± cos i D% + Di 



(31) 



where \i = is the relative proper motion of the lens with respect to the source star (as seen from the 
observer) and v± = \v± \ is the relative transverse linear velocity of the lens. 

If the caustic line is at x = 0, and the center of the star is at x = a < 0, then the amplification profile 
of the partial images (a > — 1) is given as follows. 

. . Jf n A x dx 1 f 1 (l-x 2 )^ , 

9n (a) = J Jn x = / V ' dx (32) 

J /„ dx 2b n+ i J a \/x - a 

For the two full extra images of the star when it is completely inside the caustic curve (o < —1), the 
integration range should be extended to [—1,1]. The functions g n (a) in figure |3] and ^| show some of the 
characteristic effects of the luminosity profiles at a line caustic crossing: 



1 . The lightcurve becomes practically insensitive to the luminosity profile when the star is about two 
stellar radii away from the caustic line (a = —2). 

2. In the immediate neighborhood of a = — 1, the lightcurve rises rapidly and the curvature of the light 
curve changes from concave upward to convex. Thus, given a well-sampled lightcurve, the caustic 
entry "time" a = — 1, where the limb of the star touches the caustic line producing two images joined 
at a critical point, can be read off from the lightcurve with a relatively small error. Past a = — 1, 
the images become partial, and the loss of the flux due to the partial imaging causes the slope of the 
lightcurve to decrease leading to peak turn-over and rapid decline of the lightcurve. 

3. The peaks are relatively broad, and the practically universal intersection points of the falling curves at 
a = a/ w 0.28 offer a natural distinction between the peak (a < a/) and the tail (a/ < a < 1). With 
increasing n (the surface luminosity is more concentrated toward the center of the star with larger n), 
the rising curve of the peak becomes less steeper and the falling curve becomes more steeper. (The 
slope of g diverges at a = — 1.) 
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4. The peak crossing time where the lightcurve has the maximum amplification and the peak 
amplification value are sensitive to the luminosity profile. The peak crossing times are 
a = -0.656, - 0.565, - 0.503, and -0.423 for n = 0, 0.5, 1 and 2 respectively. (The peak 
crossing time of a point source is a = 0.) The varying peak crossing times look a bit more impressive 
in unnormalized lightcurves ( Rhie 1999| ). 



5. The tails (a/ < a < 1) are dominated by the "linear behavior" except very near the end point (a = 1). 
Thus, given a well-sample lightcurve, the caustic exit "time" a = 1, where the limb of the star 
completely exits the caustic and the two extra images disappear, can be more or less read off from the 
lightcurve. The slope of the lightcurve can change abruptly (n = 0) or smoothly (n = 0.5,1,2)) at 
the end-point (see below). Near the end-point, the observed lightcurve is dominated by three normal 
images whose magnification behavior outside the caustic curve crucially depends on the proximity to 
a nearby cusp. (A cusp dominates the magnification behavior of the source plane around it more like 
the point caustic of a single and may be considered as a point caustic with directionality.) If the flux 
of the normal images is about 10% of the peak flux, the shape function may not be representative of 
the behavior of the observed lightcurve when a > 0.8. In the case of MACHO-98-SMC-1, the flux 
of the normal images was about 7% of the peak flux for the second caustic crossing and about 10% 
for the first caustic crossing ( |Rhie et al. 1999| ). In the event MACHO-97-BLG-41, it was about 16% 



'o 



for the central caustic crossing and about 30% for the planetary caustic crossing (Bennett et al 1999 



Albrow et al. 1999Q . 

6. The caustic crossing time, a = 0, marks a most featureless part of the lightcurve. 
In order to examine the the end-point behavior of the shape function g n {a), we rewrite the expression 



in ( |32[ ) with £ = 1 — x. 

-t , ri , n + l 

9n(a) = / , (33) 

2o n+ i Jo v 1 - a ~ 4 

For a « 1, the higher order term £ 2 in the numerator can be ignored, and g n (a) scales in (1 — a) as follows. 

g n {a) oc (1-a)^ (34) 
where the proportional constant is given as follows. 

2o n+ i Jo y/l-rj 

Therefore, the derivative of g n vanishes at a = l unless n = 0, which can also be seen in figure Q For a 
homogeneous model, the shape function g a has a finite slope at the end-point: g' a (l) = —V2. This has an 
interesting ramification that every lightcurve of line caustic crossing should exhibit an abrupt change of the 
slope at the end-point. That is because the luminosity profile of any star has (or is believed to have) a sharp 
edge, or, equivalently, a nonvanishing c . 

In order to gauge how well one can determine the luminosity profiles, we examine the linear models. 

c g (a) +C!ffi(a) 

giAnearW = , ; C Q + C X = I (36) 

C + Cl 
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Figure |5| shows that the shape functions for c\ = 0.5 and c\ = 0.6 are almost indistinguishable in this plot. 
The intersection points are at a = —0.76 and a = 0.28, and it will be necessary to sample the lightcurves 
in all three regions defined by the intersection points. The derivatives of the shape functions are shown in 
figure ^. The curvature changes at a = — 1 marks the entry time, and it will be relatively easy to reconstruct 
the correct time in a fitting. The linear models with c\ = 0.5 and c\ = 0.6 show relatively constant slopes 
between a = 0.3 and 0.8 or so. The opposite trends of the slopes of g and g\ make the slope of the shape 
function of a typical main sequence star with a medium value of the linear limb-darkening parameter (say, 
ci = 0.3 — 0.7) more or less constant. This "linear behavior" of the lightcurve in the tail was measured for 
the first time in the event MACHO-98-SMC-1 ( |Afonso et al. 1998| ). 

Figure shows the differences 5g Linear {a : c\) = gunear(a : ci + 5ci) - gunearifl : ci) : da = 0.1 
for ci = 0, 0.1, 0.9. In this figure, the three regions divided by the intersection points are manifest. 
The larger the value of c\, the higher the peak value at a = —0.28 and the larger amplitude of deviation. 
Thus, the smaller the value of c\ , the more demanding is the photometric accuracy and precision. In order 
to see the level of photometric accuracy we need, we normalize the difference curve 5 g Linear ( a '■ c\ ) by 
9Linear{d : c\ + 0.5(5ci) + 0.2 as an example. Figure [8] shows the normalized difference curves. In order to 
determine the linear limb-darkening parameter c\ , one needs to be able to reconstruct the shape function, 
and this requires that all three regions defined by the intersections points a = —0.76 and a = 0.28 be 
sampled fairly. The differences in tails look more discriminatory. However, the tails are more susceptible 
to the effects of the normal images and blending, and should be interpreted properly case by case. In 
an observed lightcurve, the determination of the scale factors A c and A* can be compromised if one 
overinterprets the "linear behavior" of the tails. 

In summary, we have analyzed the magnification behavior of the two extra images that appear and 
disappear in a line caustic crossing binary microlensing. We have identified the shape function of the 
lightcurve of the two extra images and examined its dependence on the luminosity profile of the lensed 
star. We also discussed the multiplicative factors that determines the overall amplification and the ones that 
scale the time axis in an observed lightcurve. We suggest that it is important to sample the lightcurve fairly 
in all three regions that are defined by the intersection points. The span of the lightcurve that is affected by 
the luminosity profile is about 3r* (a = [—2, 1]). In the case of linear models, we estimate from figure [j§ 
that the photometric accuracy of 0.3-0.8 % is required to be able to determine the linear limb-darkening 
parameter c\ with precision of 5c\ = 0.1. It is similar to the estimation made for the eclipsing binary 
method by Popper (1984). In lensing, the signal lies in magnification while it is in absorption in the case of 
eclipsing binaries. Thus, microlensing will have advantage for fainter stars. 
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2(f 

Fig. 1. — y/\d-J\ as a function of the phase angle of k = e l2v : The central caustic has four cusps, and 
\d-J\ has four zeros. Note that the total change of the phase angle along the caustic curve with four cusps is 
4-7T. (It is 2it for a triangular caustic curve and 8it for a hexagonal caustic curve. The total sum of the phase 
angle change along all the caustic loops of a binary lens is 8tt which is a topological invariant related to the 
number of point masses that is two for a binary lens.) 
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Fig. 2. — The luminosity profile functions across the diameter of the stellar disk (upper panel) and 1-d 
luminosity profile functions (lower panel) for n = 0, 0.5, 1, 2: The total flux emitted perpendicular to the 
disk (or toward an observer) has been normalized to it. In a line caustic crossing microlensing, one can 
effectively scan the 1-d luminosity profile in time. 
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Fig. 3. — Light curves of the two extra images of stars with luminosity profiles given by f n : n = 0, 0.5, 1, 2. 
a is the position of the center of the stellar disk with respect to the caustic line in units of the radius of the 
star r a st. The caustic line touches the limb of the star from inside a caustic at a = — 1 and leaves the caustic 
at a = 1. The luminosity profiles are indistinguishable when the star inside the caustic line is about one 
diameter away from the caustic (a = —2). The peak crossing times are sensitive to the luminosity profiles. 
The peak crossing times are at a = —0.656, — 0.565, — 0.503, and —0.423 for n = 0, 0.5, 1 and 2. 
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Fig. 4. — The derivatives g' n (a): at a = — 1, the slope of the shape functions g n (a) change singularly. At 
the end-point (a = 1), the slope is discontinuous for n = 0, and continuous for the others. 
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Fig. 5. — Light curves of the linear limb-darkening models: g Q and g\ determine the intersection points 
a = —0.76 and 0.28, and the shape functions of linear models distributed in between g and g\. The typical 
linear limb-darkening parameters c\ = 0.5 and c\ = 0.6 are indistinguishable in this plot. 
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Fig. 7. — Difference lightcurves for linear models with c\ values that differ by 5c% = 0.1: <%Lj near (a : 
c\ + 5c\) — 6gLinear{a ■ c\) : 5c\ = 0.1 for ci = 0.0, ...,0.9. The difference lightcurve with the largest 
amplitude is for c\ = 0.9. 
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Fig. 8. — Normalized difference lightcurves for linear models with c\ values that differ by 8c\ = 0.1: the 
difference lightcurves in figure 7 have been divided by g(a : c\ + 0.55ci) + 0.2. The baseline value 0.2 
we chose here as an example amounts to ps 12% of the peak amplification. In order to determine the linear 
limb-darkening parameter c\ with 8c\ = 0.1, photometric accuracy of 0.3-0.8 % is needed. 



